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We study cosmology in a five-dimensional brane- world with a stabilizing eflFective potential for the 
radion and matter localized on the two branes. We consider the corrections induced by the Gauss- 
Bonnet contribution to the total action performing an expansion around the two possible static 
solutions up to second order in the ratio between brane matter energy density and brane tensions. 
The Friedmann and acceleration equations on the visible brane are obtained and discussed. 

PACS numbers: 



I. INTRODUCTION 

Following the approach used in Q we analyze the 
cosmoloeical implications of the GB action in a RS 
setup (II H S H li)with radion stabilization through 
the addition of an effective term for the radion field 
to the total action. The Einstein equations are solved, 
assuming the adiabatic evolution of the fluids on the 
branes, up to the second order in p/V where p is the 
energy density of the fluids and V the tension of the 
branes. Finally the corrections to the Friedmann and 
acceleration equations are analyzed and compared to 
the case in which the GB term is switched off. 



II. EINSTEIN EQUATIONS AND STATIC 
SOLUTIONS 



Let us start with the following action for the five- 
dimensional bulk dynamics 



Sbulk — f <^X\/~-g {C,EH + l^^GB 



where 



Ceh = R + 2k^{K-U) 



i?2 - AR^^Rab 



R^^'"^ Rabcd 



(1) 



(2) 



(3) 



with U the effective contribution which stabilizes the size 
of the extra-dimension (i^J). The boundary conditions 
for the above action are fixed by the brane contributions 
to the stress-energy tensor; in order to recover RS I when 
the GB contribution is switched off and no matter is 
present on the boundaries we consider two branes with 
different tensions located at the fifth dimension orbifold 
fixed points. If we choose the following ansatz for the 
metric 



we can express the stabilizing contribution as 



U = M{h{y,t)-rf 



(5) 



where r is the expectation value of the radion field and 
the boundary terms are given by 



diag(Vi + Pi, Vi - pi,Vi - pi,Vi - pi) (6) 



where i = p,n, yi — 0,1/2 are the positions of the branes 
and Vi's are the brane tensions. 

By varying the total action one obtains the Einstein equa- 
tions 



Gab + aHAB = fc' ^9ab + Tab + Y.^T,)ab (7) 



where Hab is the second order Lovelock tensor and 
= -diag(C/, U, U,U,U + bdU/db) is the radion po- 
tential contribution. The equations Q, given the ansatz 
form a system of four differential equations. Due to 
the fact Gab and Hab satisfy the Bianchi and Bach- 
Lanczos identities, respectively only three equations are 
independent. The bulk dynamics is thus determined by 
the system 



Gtt + aHu = fc2 (a<?; 
b' = {r-h) (2: + 3^ + 2^ 



(8) 



Gyy + OiHyy = t 



yy ' ^yy 



a\y,t)dx'Ax' + b\y,t)dy^ (4) 



(see |i[T3) where a prime denotes a derivative with re- 
spect to y. The boundary equations give the following 
junction conditions on the two branes: 
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52 |l + 2«- 




(9) 



r 



where a dot denotes a derivative with respect to the uni- 
versal time t. Note that the eqs. (|8I9II reduce to the 
standard case (without GB contribution) in the Umit 
a — > 0. Before investigating cosmology, in the setup just 
described, for pi ^ 0, we should consider the case pi = 
in order to find the static solutions to perturb about. It 
is easy to verify that a warped metric still satisfies eqs.(|Sl- 
O. In fact, if one makes the ansatz 

a{y,t) = n{y,t) = exp[-mry]; b{y,t)^r (10) 

the Einstein equations and junction conditions are veri- 
fied provided 



m = ± 



with Ifc^aA < 1 and 




6m / 2 „ 



(11) 



(12) 



that a tuning of brane tension is required only if one looks 
for static solutions: different tunings could be, as usual, 
treated as a perturbation pi — AV^ and pi = —AT/; and 
generate an expanding (contracting) phase which can be 
studied by means of the formalism we introduce in the 
next section. 



III. THE PERTURBED SOLUTIONS 

We now investigate the cosmological evolution in the 
brane-world with GB contribution. The calculations are 
based on the perturbative approach {pi/Vi ^ 1) dis- 
cussed in details in ,1] (starting from a slightly different 
set of equations). The starting point is a static solution 
of the form (|10|l where 13j m — \moid\ or m — \mnew\- 
When some kind of matter is added to the branes the 
system is detuned and evolves in time. The solution be- 
comes 



n{y,t) = exp{-mr\y\) [l + 6fn{y,t)] 



(15a) 



Apart from the two expected solutions which reduce to 
the static RS when a ^ 0, namely 



a{y,t)^ah{t) expi-mr\y\) [1 + Sfaiy,t)] (15b) 



mold = ± 



two additional solutions are obtained (see 11]) which are 
interestingly less sensitive to the bulk cosmological con- 
stant, since in the limit a — > the warp factor becomes 
independent of the bulk content 



= ± 



\ 




Furthermore, due to the GB terms, warped solutions are 
still present when the bulk is filled with a positive cos- 
mological constant as (|11|) is still real in that case. Note 



b{y,t)^r + SMy,t) 



(15c) 



and the differential equations © can be written in terms 
of Hh = fl/i/a/i, 6fp (with f3 — a,n,b) and their deriva- 
tives. Moreover we note that when such a detuning is 
small, one can rely on the perturbative ansatz 



Hi 



f0.liy)Pp + f0.niy) Pn 



(1), 



(16a) 



ff^iiy) pI + ffliy) pi + fflaiy) pp Pn 



y^h.p Pp ^ "'h,n Z^" ) 

2 2 
h,p Pp 



^hln Pp Pn , 



(16b) 
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and expand these equations up to second order in e. Note 
that the time dependence in the approximate solution 
above is encoded in pi{t) which evolves, as usual, sat- 
isfying the constraint of the continuity equation for a 
4-dimensional fluid with equation of state pi — WiPi that 
is Pi = -5[d{yi,t)/a{yi,t)]{l + Wi)pi. One can finally, or- 
der by order in e, solve eqs. 0, which now contain just 
derivatives with respect to y, and determine the integra- 
tion constants and the unknowns parameters in (jl6b|) by 
making use of the boundary conditions © • Iterating the 
procedure described above up to second order is almost 
straightforward (see m for details). Due to the fact that 
the Pi's evolve independently, one needs to solve 15 dif- 
ferential equations for the coefficients /i^?, fj^}. Once 
such bulk equations are solved, one is left with 25 co- 
efficient to be fixed (10 of which derive from first or- 
der calculations) and the junctions conditions still to be 
imposed. It is possible to divide these coefficients into 
two categories: the ones which are related to the gauge 



freedom of the metric (they have to do with the defini- 
tion of the time and of the three-dimensional scale factor 
a{t)), and the ones that are related to the boundary dy- 
namics. The junction conditions form a system of 15 
independent equations which determine just the dynam- 
ical coefficients. Five gauge coefficients can be fixed by 
arbitrarily defining a time evolution parameter and five 
(those related to the scale factor) may remain arbitrary 
since they are not present in the expressions of physical 
observables. 

The Friedmann and the acceleration equations on the 
negative tension brane (the negative tension brane, in 
RS I setup with m > 0, is supposed to be the 4- 
dimensional space-time manifold in which we live) can 
thus be obtained by fixing n{l/2,t) = 1 and expressing 
= (d(l/2,i)/a(l/2,t))2 and &/a = a(l/2, t)/a(l/2, i) 
as functions of the time measured on the same brane. 
This leads to (we let e — > 1 at the end of the calcula- 
tions) 



a _ Pm [{3wp + l)e2"Vp + (3w„ + l)p„] 27k^m^ (e"" + 1) (w„ + 1) (w„ + |) (w„ - |) ^ 
a ^ die-^^ -l){2m^a + l) 32Mr'^ (e""^ - if {2m^a + 1) ^" 

\ , ^, , fc^m^a (e™' + 3) (3w.p + 1) (3i«„ + 1) m^k^ (l - Am'^a^) 

+ \ — rr(3w„ -I- l)(3w„ -I- 1) H ——, ^ r 

\l44^ " ' 72(e'"'--l) 96r2M (e""- - 1)' 

X [e2"'^(3w;„ - l)(18w^ -I- 9wpW„ -f Ihwp + 39w„ + 23) + 9e"'^ {<owl + &wl + 3w^w„ 

+iWpWj^ -\- llWp + llWn -I- 2WpWn + Wp + W„ — 4] + MWp + 27WpWn + 99Wp — ISWpWn 

e"^^ppp,, [ m2fc^(3^;p-l)(l-4mW)(e""'-H) 
-|-24wp — 21w„ — 13J ^ 



{2m^a-l){2m^a + lf \ 96r2M (e™'' - 1)^ 



r mr/n^ 2 r. noN n rl ^^^/n 1^2 A^^m^a (1 + 3e"'') (3u.p -f 1)^ 

X [e^^l27wl + Mwp + 23) - 9.«p - 5] + — (3t.p + 1)^ (e- - 1) 



2mr 2 
Hp 



{2rn?a-^l){2m^a + lY 



(17) 



khn , 3fc2^2 (1 ^ gm.) ^ ^) ^^^^ _ 1^ ^ f fcV^g (e""- + 3) 

3 (e"'- - 1) {2m^a + 1) ^ ^"^ 32r2M (e"- - 1)' {2m^a + 1) ^" I 18 (e"'^ - 1) 

k^ m^/c^ (1 - 4m4a2) (e""- + 1) r ? . s mr, ■? ^^ 1 

+ - + 48,2M(e^._\)2 [18-^ - ~ 1) ~ - 1 + 2e"'^(9-^ + " 4)] | 



(2m2a - 1) (2TO2a + 1)^ \ 96r2M (e™'' - 1) 

k^ k^m^aiie^'-' + \)\ e^"'"- pi 



36 I8(e™'--1) J (2m2a- l)(2m2a + l)'' 

(18) 



The above expressions are quite involved as they contain the ones already found in |llll2| for a — s- and m — moid- 
contributes from the fluids on both branes and reduce to 
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Note that, due to the choice of the time parameter, the 
coefficient of p^-^ vanishes when p„ behaves as vacuum 
energy or radiation. As a consequence, when pp is negh- 
gible, the second order terms, which would be otherwise 
responsible for a deviation from the standard cosmolog- 
ical evolution, vanish for Wn = —1 or Wn — 1/3. This 
important feature of brane-world scenarios with radion 
stabilization was already present in the case a = and is 
furthermore conserved when GB contribution is present. 
After some algebra one can partially absorb the GB cou- 
pling in lfTFI) - ((T7|l by redefining fc^ = k"^ / {im^ a+1) which 
means that GB corrections can be observed only through 
the indirect contribution of positive tension brane mat- 
ter. Furthermore note that a is always multiplied by 
m? and consequently the GB contributions vanish in the 
limit a — > only if we consider the case m = moid- In 



the other case m 



lim a 



one has 



= 1 



(19) 



and the evolution obeys a dynamics which is modified 
with respect to the standard RS case. In this case, keep- 
ing just the leading contributions and letting Wp — > — 1 
(small brane detuning), one is left with 



fc2 r 
9 



mpp - (3w„ -I- 4) (3w„ - 1) 



, m^PpPn 

8r^M 



(20) 



a 
a 
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2mpp + {3wn — 1) {3wn + 4) 



X (3w„ -f 1) 



m^PpPn 
8r^M 



m^Pp 



(21) 



where pp = e^^pp and one should keep m • pp and m • p„ 
small. Apart from the contribution of pp to the effec- 
tive cosmological constant, the term proportional to p^ 



is negligible and one observes the usual contribution to 
the expansion rate at Wn = — 1 or Wn = 0. 
Finally we note that, when M — > oo, some of the sec- 
ond order terms vanish, as they did in the standard case 
a = 0. In fact these terms derive from the radion dy- 
namics; such dynamics is sensible to the state equation 
of the fluids on the branes and generates a complicated 
Wi dependence in (|18I17|) . In the limit M — > oo when the 
radion is fixed to the minimum of the stabilizing potential 
r it becomes trivial and these terms vanish. 



IV. CONCLUSIONS 

We have examined the four-dimensional cosmologi- 
cal equations deriving from the Einstein equations for a 
brane-world with a stabilizing potential in the presence 
of a Gauss-Bonnet term in the action. We found solu- 
tions performing an expansion up to the second order in 
p/V in order to examine the cosmological behavior on 
the negative tension brane. The formalism can be easily 
extended to describe the positive tension brane as well. 
Due to the Gauss-Bonnet extra terms, the system admits 
two different static solutions: one behaves as the usual 
RS when the Gauss-Bonnet coupling a goes to zero, while 
the second one has a warping factor independent of the 
bulk content in the limit a — > . 

At first order the deviations from standard (without GB 
term) equations can be reabsorbed with a redefinition 
of the 4-dimensional Newton constant. The same holds 
when one considers second order equations with a van- 
ishing Pp. On the other hand, if one considers the con- 
tributions due to matter perturbations on the positive 
tension brane some deviations appear. An interesting 
feature emerging is that in the limit a —^ such devia- 
tions are not swept away when perturbing the solution 
that does not reduce to the usual RS. 

We would like to thank R. Casadio for the helpful com- 
ments. 
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